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Charm production from the quark-gluon plasma created in the midrapidity of central heavy ion collisions at 
the Large Hadron Collider (LHC) is studied in the next-to-leading order in QCD. Using a schematic longitudi- 
nally boost-invariant and transversally expanding fire-cylinder model, we find that charm production could be 
appreciably enhanced at LHC as a result of the high temperature that is expected to be reached in the produced 
quark-gluon plasma. Sensitivities of our results to the number of charm quark pairs produced from initial hard 
scattering, the initial thermalization time and temperature of the quark-gluon plasma, and the charm quark mass 
are also studied. 
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I. INTRODUCTION 

In heavy-ion collisions at relativistic energies, hadrons 
composed of confined quarks and gluons are expected to dis- 
solve into their constituents and form an extended volume of 
quark-gluon plasma (QGP). Experiments at the Relativistic 
Heavy Ion Collider (RHIC) have indeed shown that the re- 
sults are consistent with the formation of a strongly interact- 
ing quark;^luon plasma during the initial stage of the colli- 
sions mUlllHl- Many observables have been proposed to 
probe the properties of the quark-gluon plasma. Among them 
is the production of particles consisting of heavy charm and 
bottom quarks. In particular, it was suggested that the produc- 
tion of charmonium J/\\t in relativistic heavy ion collisions 
might be suppressed as a result of its dissociation in the pro- 
duced quark-gluon plasma due to the Debye screening 
Although recent studies based on the lattice QCD have in- 
dicated that the J/\\t can survive in the quark-gluon plasma 
at temperatures up to about twice the deconfinement temper- 
ature 161 170, J/\\i production in relativistic heavy ion colli- 
sions may still be suppressed in heavy ion collisions at RHIC 
lllllilll^l. However, if the initial produced charm quark pairs 
are large such as in heavy ion collisions at LHC, pro- 
duction could be enhanced through regeneration from charm 
and anticharm quarks in the quark-gluon plasma ifTlll as well 
as through statistical production from lfl2|] or coalescence of 
|[l3ll charm and anticharm quarks during hadronization of the 
quark-gluon plasma. These mechanisms for J/\\t production 
would become even more important if there are other sources 
for charm quark production in relativistic heavy ion collisions. 
To use charmonium and bottomonium production in relativis- 
tic heavy ion collisions as a diagnostic tool for the properties 
of produced quark-gluon plasma ||9l[l^[Tll[l3l, it is thus es- 
sential to understand the production mechanism of charm and 
bottom quarks in these collisions. 

In heavy ion collisions, there are generally four different 
contributions to charm production: direct and prethermal pro- 
duction from partonic interactions as well as thermal produc- 
tion from partonic and hadronic interactions. The direct pro- 
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duction is from initial hard scatterings of partons in the nucle- 
ons from the colliding nuclei, which happens at a time scale 
of about 1/OT( ~ 0.15 fm/c with being the charm quark 
mass. It contributes about 3 charm quark pairs in one unit 
of midrapidity in AuH-Au collisions at ^/snn = 200 GeV at 
RHIC, and this number is increased to about 20 in Pb-nPb 
collisions at = 5.5 TeV at LHC fillill. These num- 

bers have, however, substantial uncertainty, particularly in 
extrapolating to LHC energy from that at RHIC lfT2l [Till , 
where measurements from PHENIX for pH-p collisions jlql 
and STAR for dn-Au collisions [|l7i1 disagree by about a factor 
of 2. The prethermal production is from produced partonic 
matter before it reaches thermal equilibrium. Based on con- 
sideration of minijet partons, it has been shown in Ref. lfisll 
that this contribution is unimportant compared to that from 
direct production. For thermal production, it can be from the 
thermalized partonic matter or quark-gluon plasma as well as 
from the hadronic matter formed after the hadronization of 
the quark-gluon plasma. For thermal production of charm 
quark pairs from the quark-gluon plasma, studies so far have 
been based on the lowest-order QCD process of gluon-gluon 
fusion. The results show that its contribution depends sen- 
sitively on the initial temperature of produced quark-gluon 
plasma lfl9[ l20l l2ll l22i l23ll . and could be important if the 
initial temperature is high. In hadronic matter, charmed 
hadron can be produced from reactions such as 7tA^ DAc 
iH, pA? DAc iH, and NN NAcD il. According to 
Ref. luTIl based on the Hadron-String Dynamics (HSD) with- 
out a partonic stage, secondary meson-baryon reactions in the 
hadronic matter gives rise to approximately 11 % enhancement 
of charm pair production in heavy ion collisions at RHIC en- 
ergies. This contribution is, however, expected to be much 
smaller if one takes into account the formation of the quark- 
gluon plasma in the collisions, which would lead to a final 
hadronic matter with significantly lower density and tempera- 
ture than in the HSD model. 

For AuH-Au collisions at y^s^N — 200 GeV at RHIC, charm 
production is dominated by direct production as the contri- 
bution from thermal production is negligible due to the rel- 
atively low initial temperature (~ 350 MeV) of produced 
quark-gluon plasma. Since the thermal production rate in- 
creases exponentially with temperature as we will show later, 
thermal production of charm quark pairs in Pb-nPb collisions 
at ^/snn — 5.5 TeV at LHC, in which a quark-gluon plasma 
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with much higher initial temperature and density than those at 
RHIC is expected to be created, will be greatly enhanced. On 
the other hand, charm production from initial hard scattering 
of colliding nucleons increases only logarithmically with the 
collision energy. In most studies on the production of charm 
particles at LHC, it has been assumed, however, that thermal 
production can be neglected as well, and only the direct pro- 
duction from initial hard scattering has been considered. In 
the present study, we will consider thermal charm production 
at LHC and compare its contribution to that due to direct pro- 
duction to see whether this assumption can be justified or un- 
der what conditions it can be considered as a good approxi- 
mation. 

Charm production in a quark-gluon plas ma has been previ- 
ously studied in the lowest-order QCD |[l9il20i l2ll l22ll28il29ll . 
In our study, we will consider thermal charm production in the 
next-to-leading order and investigate its sensitivity to the ini- 
tial conditions of the produced quark-gluon plasma at LHC. 
We find that based on a reasonable estimation of the initial 
conditions, thermal charm production could play an important 
role in determining the total number of charm pairs produced 
in relativistic heavy ion collisions. 

This paper is organized as follows. In the next section, 
we discuss the processes and corresponding cross sections for 
thermal charm production in the leading order as well as the 
next-to-leading order in QCD. In Section III, we introduce the 
time evolution of formed quark-gluon plasma in heavy ion 
collisions at LHC and derive the rate equation for charm pro- 
duction in the quark-gluon plasma. In Section IV, we present 
the numerical results on thermal charm production at LHC 
and their dependence on the initial conditions and other pa- 
rameters in the model. Finally, a summary and conclusions 
are given in Section V. 



n. CHARM PRODUCTION FROM QUARK AND GLUON 
INTERACTIONS 

A. charm production cross sections 

Previous studies of thermal charm production in a quark- 
gluon plasma Jil HH HI [H were all based on the fol- 
lowing leading-order processes in QCD: 

q + q c + c, (1) 
g + g c + c. (2) 

In the present study, we include also the next-leading-order 
processes 

q + q c + c + g, (3) 
g + g c + c + g, (4) 

and the interferences between the leading-order processes 
with their virtual corrections due to vertex corrections and self 
energy insertions as their contributions are of same order in 
the QCD coupling as the next-leading-order processes. We 
neglect, however, next-leading-order processes involving the 



gluon-quark interactions 

g + q c + c + q, 

g + q ^ c + c + q, (5) 

as their contributions are less important compared to other 
processes, especially the gluon-gluon fusion process of 
Eq.lHl). This is largely due to the smaller color factor from 
the gluon-quark interaction vertex as compared to that from 
the gluon-gluon vertex, and the fact that the charm pair can be 
produced from two incoming gluon lines in gluon-gluon fu- 
sion process, whereas in the gluon-quark process it can only 
come from a single gluon line ll30ll . 
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FIG. 1 : (Color online) Charm quark pair production cross sections 
as functions of center-of-mass energy from quark-antiquark annihi- 
lation (left panel) and gluon-gluon fusion (right panel) in both the 
leading (dashed line) and the next-to-leading (solid line) order. The 
charm quark mass is taken to be nic =1.3 GeV, while light quarks 
and gluons are massless. 

The total cross section for charm production in the next- 
to-leading-order in QCD has been previously derived. In the 
present study, we use the results given in Refs. Isollsill . where 
the soft divergences in the next-leading-order 2 — > 3 processes 
are eliminated by corresponding divergences in virtual correc- 
tions to the leading-order 2^2 processes, while the ultravi- 
olet divergences in virtual corrections are handled by renor- 
malization in the modified minimal-subtraction {MS) scheme 
ll32[ I33I1 . Choosing the QCD running coupling constant as 
af(mc) and the renormalization scale as the charm quark 
mass, we have calculated the cross sections for the production 
of charm pairs with charm quark mass nic = 1.3 GeV from 
massless light quarks and gluons, and they are shown in Fig.[T] 
as functions center-of-mass energy for the quark-antiquark 
annihilation (left panel) and gluon-gluon fusion (right panel) 
processes. The dashed and solid lines correspond to results for 
the leading order and the next-to-leading order, respectively. It 
is seen that for charm production from quark-antiquark anni- 
hilation, the next-to-leading order gives a larger cross section 
at low energies but a smaller one at high energies. The cross 
section for charm production from gluon-gluon fusion is, on 
the other hand, significantly larger in the next-to-leading order 
than in the leading order at all energies. 
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B. thermal averaged charm production cross sections 

In the kinetic model to be used in the next section for study- 
ing charm quark production and annihilation in a quark-gluon 
plasma, thermal averaged cross sections are needed. In terms 
of the thermal distribution functions /, (p) of quarks and glu- 
ons in the quark-gluon plasma and the relative velocity v„/, of 
two initial interacting partons a and b, the thermal averaged 
cross section Oat^ai for the reaction ab — > cd is given by ll34ll 

"^^""^ J d^Pad^Pbfa{Pir)MPb) 

= [4alK2{aa)alK2{ah)]-' 

X / dz[z^ - (a„ + abf] [z^ - (a„ - abf] 

J ZD 

X K,{z)o{s=z^T^), (6) 

with a,- = nii/T, zo = max(aa + OLh,a.c + cxj), and Ki being 
the modified Bessel function. 



III. CHARM PRODUCTION IN HEAVY ION COLLISIONS 
ATLHC 

Using above thermal averaged charm production cross sec- 
tions, we study in this section the time evolution of the abun- 
dance of charm quark pairs in heavy ion collisions at LHC 
using a kinetic model based on the rate equation that takes 
into account both production and annihilation of charm quarks 
in the produced quark-gluon plasma. For the dynamics of 
the quark-gluon plasma, we describe it by a schematic hy- 
drodynamic model and assume that both quarks and gluons 
are in thermal and chemical equilibrium during the evolution. 
We further assume that all produced charm quarks including 
those produced directly from the initial hard collisions are also 
in thermal equilibrium, although not in chemical equilibrium. 
The latter is consistent with observed large elliptic flow of the 
electrons from charmed meson decays in heavy ion collisions 
at RHIC IjtI [ssll . which requires that charm quarks interact 
strongly in the quark -gluon p lasma and thus are hkely to reach 
thermal equilibrium ll39ll4otl4lll . 
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FIG. 2: (Color online) Thermal averaged cross sections for charm 
pair production from quark-antiquark annihilation (left panel) and 
gluon-gluon fusion in a quark-gluon plasma as functions of temper- 
ature. The charm quark mass is taken to be nic =1.3 GeV while 
quarks and gluons are taken to have thermal masses given by Eq. (|7). 

In the quark-gluon plasma, quarks and gluons acquire ther- 
mal masses. As an exploratory study, we include this effect in 
their distribution functions but not in the calculation of their 
scattering cross sections. For t he therm al masses of quarks 
and gluons, they are taken to be | 



--gT/\/6 and m^= gT j^fl, 



(7) 



where g is the QCD coupling constant and is taken to have 
the value g ~ ■\/47tas(27tr). With the charm quark mass 
rtic = 1.3 GeV, the thermal averaged cross section for charm 
production in a quark-gluon plasma is shown in Fig. |2] as 
functions of the temperature of the quark-gluon plasma for 
quark-antiquark annihilation (left panel) and gluon-gluon fu- 
sion (right panel). For both reactions, thermal averaged cross 
sections are significantly larger in the next-to-leading order 
(sohd fine) than in the leading order (dashed line). 



A. the rate equation 

The time evolution of the number density of charm quark 
pairs «tr in an expanding qua rk-g luon plasma can be de- 
scribed by the rate equation 11421 14311 : 



qq^ccg 



ccg^qq - 



gg^cc - 



cc^gg ' 



gg^ccg 



i:cg^gg ■ 



(8) 



In the above, = y(1,v) is the four velocity of a fluid ele- 
ment in the quark-gluon plasma with velocity v and the cor- 
responding Lorentz factor y, and the terms on the left hand 
side of above equation are the charm pair production and an- 
nihilation rates. To the next-to-leading order, the charm pair 
production rate is given by 



^qq^cc 

Rgg^CC 

Rqq^ccg 



= {Oqq^ccV)ngng, 



1 



-{Og^^^ccgv)ng, 



(9) 



where rig, rig and rig denote the gluon, quark, and antiquark 
densities in the quark-gluon plasma, respectively, and they 
are taken to have their equilibrium values. The leading-order 
cross sections Oqq^cc and O^g^cc in above equation are com- 
puted from the processes in Eq.([T]i, while the next-leading- 
order cross sections Oq^^ccg and Ogg^ccg include both the 
processes in Eq.© and the virtual corrections to the leading- 
order processes. 



In Fig. [3] we show the thermal charm production rate as a 
function of temperature with massive quarks and gluons in a 
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FIG. 3: (Color online) Thermal charm quark production rate as a 
function of temperature for massive partons and charm quark mass 
»ic = 1.3 GeV. The inset gives the ratio of the charm production rate 
in the next-to-leading order to that in the leading order. 



thermally and chemically equilibrated quark-gluon plasma. It 
is seen that the charm production rate increases almost expo- 
nentially with increasing temperature. With the initial temper- 
ature To ~ 350 MeV in central heavy ion collisions at RHIC, 
the thermal charm production rate is more than two orders of 
magnitude smaller than that at the temperature Tq ~ 700 MeV 
expected to be reached in central heavy ion collisions at LHC. 
This thus justifies the neglect of thermal charm production at 
RHIC. The inset in Fig. |3] gives the ratio of the charm pro- 
duction rate in the next-to-leading order to that in the leading 
order, which is seen to vary from ^ 4.5 at low temperatures to 
1.8 at high temperatures. 

Since the charm pair annihilation rate and production rate 
become equal when the charm quark abundance reaches 
chemical equilibrium, i.e. 
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where n^^ is the equilibrium density while R^^ is the rate eval- 
uated with the equilibrium densities, the charm pair annihila- 
tion rates in the leading order and the next-leading order can 
thus be written as 

Rcc^qq = {<3qq^ccV){n'^ f {^^^^ , 
1 , 



Rcc 



Rtxa 



-gg = :^{(ygg^ccV)itf^) ^ ^eq 



^2 I '2£i 



g^gg 



= ^{<ygg- 



^c,gy)in^)H^ 



(11) 



In the above, we have assumed that the density of quarks is 
the same as that of antiquarks, i.e., n^'^ = n^^. 

With above equations and taking light quarks and gluons 
to be in thermal and chemical equilibrium, the rate equation 
Eq. (O can be rewritten as 

= [{{Ociq^ccV) + {(Sqq^ccgv)){n'^f 



l-( — 



(12) 



B. collision dynamics at LHC 



In central heavy-ion collisions at RHIC, the particle dis- 
tribution at mid-rapidity is approximately uniform, so the 
Bjorken boost-invariant hydrodynamic solution is a good ap- 
proximation for describing the time evolution of the hot dense 
matter formed in these collisions ll44ll . Such a model indeed 
give s a good description of the experimental data at RHIC 
14511. Also, recent studies based on the AdS/CFT correspon- 
dence have demonstrated that at late proper times the system 
produced in relativistic heavy-ion collisions in the strongly- 
coupled regime described by gauge-gravity duality exhibits 
an energy density scaling with characteristics similar to the 
Bjorken hydrodynamic solution (H^H^j. We thus expect that 
the longitudinal dynamics in heavy ion collisions at LHC is 
also boost invariant. To include the effect of transverse flow, 
we further include an accelerated transverse expansion, result- 
ing in a cylindrical symmetry in the geometry of the collision. 
In terms of the cylindrical coordinates r and (p as well as the 
proper time T and the space-time rapidity r| defined by ll42ll43ll 



2 t~z 



(13) 



we then have = u'^ = 0. If we further assume a uniform 
density distribution in the transverse plane and take the aver- 
age over the radial coordinate, the left hand side of the rate 
equation Eq. (fT2l) can then be expressed as; 



1 



{zR^{z)n,,{in), 



(14) 



where R{x) denotes the transverse radius of the system, and 
{u'^) is the averaged x component of the four velocity defined 
as 



drru^{r) . 



(15) 



At mid-rapidity, the four velocity of a fluid element can 
be described by two independent boosts in longitudinal and 
radial directions ll48ll . with the radial flow velocity pr given 
by 



= y,- 



1 



(16) 
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Assuming the usual ansatz for the radial expansion velocity 
||4l|4lS,|4il 



we have 



- dy , ■ 

y/l-{dR/dT:)2y 



(17) 



(18) 



With the time evolution of the transverse radius of the fire- 
cylinder taken to be llsoll 



R{z)=Ro + a{z-Zof/2, 



(19) 



where Rq and Tq are the initial radius and proper time of the 
fire-cylinder, respectively, and a denotes the transverse accel- 
eration, the volume of the fireball at proper time T is then given 
by 



y(T) =7t/?2(T)T 



(20) 



In the following, we choose the initial transverse radius of the 
fire-cylinder to be the radius of the colliding Pb nucleus, i.e., 
/?o = 7.0 fm, and a transverse acceleration a = 0.1 c^/fm in 
order to obtain a reasonable final transverse flow velocity. 



C. time evolution of the temperature of the quark-gluon 
plasma 

To estimate the initial temperature of produced quark-gluon 
plasma, we use two different models: A Multi-phase Trans- 
port (AMPT) model pill and the model based on the Color 
Glass Condensate 1521 ■ Using the AMPT model, we obtain 
an initial transverse energy of about 3000 GeV per unit rapid- 
ity at the midrapidity in central Pb-nPb collisions at -y/sJvw = 
5.5 TeV. Most of this initial energy is concentrated within a 
transverse radius of about 4.7 fm that is smaller than the radius 
of the colliding nuclei as a result of their surface diffuseness. 
If we take an earlier initial proper time of Tq = 0.2 fm/c than 
that at RHIC, the initial energy density of produced quark- 
gluon plasma is then about 



dEr / dy 
kRItq 



3000 



n X 4.72 X 0.2 



226GeV/fm^ (21) 



Using the relation 



2x4 



3771^ r 

30 



(r/160MeV)'*GeV/fm^ (22) 



for the energy density of a quark-gluon plasma with mass- 
less quarks and gluons, the initial temperature of the quark- 
gluon plasma formed at LHC is thus 7b ~ 620 MeV. This 
number turns out to be very close to that predicted by the 
Color Glass Condensate model. As shown in Ref.Js^, the 
glasma formed after Pb-nPb collisions at ^/snn =5.5 TeV 
has an energy density of about £ ~ 700 GeV/fm^ at a proper 



time T 0.07 fm/c. The glasma then evolves into a ther- 
malized quark-gluon plasma at a proper time Tq- Assuming 
that the energy density decreases inversely with the proper 
time during this stage, its value at Tq ~ 0.2 fm/c is then 
£0 ~ 245 GeV/fm , corresponding to an initial temperature 
Tq w 633 MeV. Since the result from the Color Glass Con- 
densate model depends on the fourth power of the saturation 
momentum, the predicted initial temperature thus has a large 
uncertainty. The same is true for the AMPT model as the pre- 
dicted total transverse energy is sensitive to the parton struc- 
ture function of the nucleon, the initial nuclear shadowing, and 
the saturation of produced gluons. Therefore, we will consider 
in the following also the scenario with an initial energy den- 
sity which is 50% or a factor of two larger larger than above 
value, corresponding to an initial temperature of about 700 or 
750 MeV, respectively. 

We note for central AuH-Au collisions at ^^Inn ~ 200 GeV 
at RHIC, the total transverse energy from the AMPT model 
is about 1000 GeV. Taking a proper time To = 0.5 fm/c for 
the formation of an equilibrated quark-gluon plasma as re- 
quired by the hydrodynamic model to explain the observed 
hadron elliptic flows [53| 113, HH, the initial energy density 
and temperature of produced quark-gluon plasma are then 
£0 w 33 GeV/fm^ and 7b « 383 MeV, respectively. These 
values are again similar to those from the Color Glass Conden- 
sate model, which gives an energy density £ w 130 GeV/fm^ 
at To = 0.1 fm/c, leading thus to £0 ~ 26 GeV/fm and 
To « 361 MeV at Tq = 0.5 fm/c. With such a low initial 
temperature, charm production from the quark-gluon plasma 
is thus negligible according to the production rates shown in 
Fig.E] 
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FIG. 4: (Color online) Time evolution of the temperature of the 
quark-gluon plasma formed in Pb+Pb collisions at ^/sivA' = 5.5 TeV 
for different initial temperatures and same initial proper time To = 
0.2 fm/c. 

For the time evolution of the temperature of produced 
quark-gluon plasma at LHC, we determine it using the entropy 
conservation, and the results are shown in Fig. |4] for the three 
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initial temperatures 7b = 630, 700, and 750 MeV. In all cases, 
the critical temperature for the quark-gluon plasma phase tran- 
sition to the hadronic matter is taken to be 7c = 170 MeV. It 
is seen that the temperature decreases very quickly with the 
proper time for all three initial temperatures, dropping below 
400 MeV after about 0.5, 1, and 1.5 fm/c for To = 630, 700, 
and 750 MeV, respectively. 



IV. RESULTS 

For central Pb+Pb collisions at ^/s^n = 5.5 TeV at LHC, 
the number of charm pairs produced in one unit rapidity from 
initial hard nucleon-nucleon collisions is about 20 at midra- 
pidity according to the next-to-leading order pQCD calcula- 
tions 11211 . Neglecting the contribution from pre-thermal pro- 
duction, we solve the rate equation with this initial number of 
charm quarks pairs and certain initial temperature and proper 
time for the produced quark-gluon plasma. We will vary the 
initial temperature and proper time to study their effects on 
the final number of charm quark pairs produced in these colli- 
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FIG. 5: (Color online) Number of charm pairs as a function of proper 
time in central Pb-l-Pb collisions at y^Smv = 5.5 TeV in the leading 
order (dashed line) and the next-to-leading order (solid line) in QCD. 

We first consider the case of an initial temperature Tq = 
700 MeV and proper time To = 0.2 fm/c. For the charm quark 
mass, we take it to be = 1 .3 GeV as that extracted from the 
experimental data in e+e^ and pp collisions. To take into ac- 
count the medium effect on quarks and gluons, we use the 
thermal masses given in Eq.(|7]i. The time evolution of the to- 
tal number of charm quark pairs obtained from these initial 
conditions in central Pb-nPb collisions at ^/snn =5.5 TeV is 
shown in Fig. |5] It is seen that including thermal production 
enhances the total number of final charm pairs as compared 
to that from initial direct production. The charm pair num- 
ber reaches the peak value at T ^ 2 fm/c and then deceases 
with the proper time. At the critical temperature 7c = 170 



MeV when the proper time is Tc ~ 6.4 fm/c, the number of 
charm pairs is about 27 in the next-to-leading order (solid 
line), which is about 30% larger than the number due to di- 
rect production from initial hard collisions. The final charm 
pair number is reduced to about 22 if only thermal production 
from leading-order contribution is included as shown by the 
dashed line. 

We note that the final charm quark pair number is much 
larger than its chemically equilibrium value at Tq = 170 MeV, 
which is about 4.7. This is not surprising as the time for charm 
pairs to reach their equilibrium number n^^ at certain tem- 
perature, given by Teq = «^.^/2/?„- with 7?„- denoting the pro- 
duction rate at that temperature, increases dramatically with 
decreasing quark-gluon plasma temperature from a value of 
a few fm/c at T = 700 MeV to a few thousands of fm/c at 
Tc = 170 MeV. It is thus not possible for charm quarks to 
reach chemical equilibrium during the finite lifetime of the 
produced quark-gluon plasma at LHC. Since the initial num- 
ber of charm quark pairs produced from hard scattering of ini- 
tial nucleons is large at LHC and the charm production rate is 
also larger during the early stage of the quark-gluon plasma 
when its temperature is high, there are more charm pairs in 
the quark-gluon plasma than the equilibrium number at the 
critical temperature. As the quark-gluon plasma expands and 
cools, the charm annihilation rate decreases, making it less 
likely to destroy the produced charm pairs and leading thus 
to an over-saturation of the charm abundance at the critical 
temperature. 
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FIG. 6: (Color online) Number of charm pairs as a function of proper 
time in central Pb+Pb collisions at ^/sjm = 5.5 TeV for different 
number of charm quark pairs produced from initial hard scattering. 

Since there is substantial uncertainty in the charm quark 
pairs produced from initial direct production lfl2i[l4ll . we have 
also studied thermal charm production using different initial 
charm quark pair numbers. Varying the initial charm pair 
number by a factor of two, we have repeated above calcu- 
lations, and the results based on the next-to-leading order in 
QCD are shown in Fig.|6] For initial numbers of charm quark 
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pairs of 10, 20 and 40, the final numbers are found to be about 
19, 27, and 45, respectively. Thermal production of charm 
quarks from the quark-gluon plasma thus becomes more im- 
portant as the initial charm pair number becomes smaller 

Returning to the case of 20 initial charm quark pairs, we 
note that in above calculations, we have used massive quarks 
and gluons given by the thermal QCD calculations. To see 
how these masses affect thermal charm production, we have 
carried out similar calculations as above with massless quarks 
and gluons but same initial conditions of To ~ 700 MeV and 
To = 0.2 fm/c as well as the charm quark mass of nic = 
1.3 GeV. We find that the resulting number of charm quark 
pairs at 7c = 170 MeV differs hardly from the massive case. 
The reason for this is that although using massless quarks and 
gluons reduces the thermal averaged charm production cross 
sections, the effect is largely compensated by the larger quark 
and gluon densities if they are massless. 



1 1 1 1 1 1 1 1 1 1 1 1 1 


/ ✓ 




/ / 




. // 










yO.2 fm/c, m^=1.3 GeV 




— T|, = 750MeV 




-- T|, = 700MeV 




.-. T|j = 630MeV 


, 1 , 


1 1 . 



1 2 3 4 5 6 
X (fm/c) 



FIG. 7: (Color online) Total number of charm pairs as a function of 
proper time in central Pb+Pb collisions at ^/sjm = 5.5 TeV for dif- 
ferent initial temperatures of the quark-gluon plasma but same initial 
proper time To = 0.2 fm/c. 

As previously mentioned, the charm production rate ex- 
hibits an exponential increase with the temperature of the 
quark-gluon plasma. The total number of charm quark pairs 
produced in heavy ion collisions thus depends on the initial 
temperature of the expanding quark-gluon plasma. In Fig. |7] 
we show the total number of charm pairs as a function of 
proper time in central Pb-nPb collisions at ^/snn =5.5 TeV 
for different initial temperatures of the quark-gluon plasma 
but same initial proper time To = 0.2 fm/c. It is seen that 
the final number of charm quark pairs at the critical tempera- 
ture Tc= 170 MeV calculated in the next-to-leading-order de- 
creases to about 22 for Tq ~ 630 MeV and increases to about 
33 for 7b = 750 MeV. 

We have also considered the effect of initial proper time on 
thermal charm production by using To = 0.5 fm/c as that in 
heavy ion collisions at RHIC. The predicted initial temper- 
ature from the Color Glass Condensate model is then Tq = 



500 MeV. Results for this initial temperature as well as those 
for the temperatures of 560, and 600 MeV, which correspond, 
respectively, to an increase of the initial energy by 50% and 
100% of that for 7b = 500 MeV, show that the final total num- 
ber of charm quark pairs is only slightly reduced compared to 
that for an initial proper time To = 0.2 fm/c, i.e., 21, 26, and 
32 for above three temperatures. 

The charm quark mass also affects the contribution of ther- 
mal charm production. With nic = 1.5 GeV and the initial 
proper time To = 0.2 fm/c, the final total number of charm 
quark pairs is reduced to about 20, 22 and 25 for the three ini- 
tial temperatures Tq = 630, 700, and 750 MeV, respectively. 



V. SUMMARY AND DISCUSSIONS 

In this paper, we have carried out the first calculation of 
thermal charm production at next-to-leading-order in QCD. 
Modeling central heavy ion collisions at LHC by a schematic 
longitudinally boost invariant and transversely expanding fire- 
cylinder of quark-gluon plasma, we have evaluated the num- 
ber of charm quark pairs produced in these collisions. With 
an initial temperature of 700 MeV for an equilibrated quark- 
gluon plasma at an initial proper time of 0.2 fm/c and a charm 
quark mass of 1 .3 GeV, we have obtained about 30% enhance- 
ment in the production of charm quarks than that produced 
directly from initial hard collisions, if the latter is taken to 
be 20 pairs at midrapidity according to the next-to-leading 
order pQCD calculations. About equal contributions are ob- 
tained from the leading order and the next-leading order pro- 
cesses. This result is, however, sensitive to the initial con- 
ditions for the produced quark-gluon plasma as well as the 
charm quark mass. The enhancement is increased to about 
80% if the initial temperature is increased to 750 MeV but re- 
duced to about 10% if the initial temperature is decreased to 
630 MeV. Delaying the proper time at which a thermalized 
quark-gluon plasma is formed does not affect, however, much 
thermal charm quark production as the effect due to decreased 
initial temperature is compensated by that from the increased 
volume of the quark-gluon plasma. Changing the charm quark 
mass has, on the other hand, a large effect on thermal charm 
quark production from the quark-gluon plasma. With a larger 
charm quark mass of nic = 1.5 GeV, thermal charm quark pro- 
duction from the quark-gluon plasma becomes unimportant 
even for an initial temperature of 7b = 700 MeV. Finally, ther- 
mal charm production from the quark-gluon plasma becomes 
more important if the number of directly produced charm pairs 
from initial hard scattering is small. These results are not only 
of interest in their own right but also useful for understanding 
charmonium production in relativistic heavy ion collisions as 
its production from the quark-gluon plasma is proportional to 
the square of the charm quark numbers. 
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